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The angular distributions of the decay products in the successive decays X^ixt) J/V'(T) 7 and 
J/V'(T) — >■ £^£~ are calculated as a function of the angular momentum composition of the decaying 
X meson and of the multipole structure of the photon radiation, using a formalism independent 
of production mechanisms and polarization frames. The polarizations of the x states produced in 
high energy collisions can be derived from the dilepton decay distributions of the daughter J/i/) or T 
mesons, with a reduced dependence on the details of the photon reconstruction or simulation. More- 
over, this method eliminates the dependence of the polarization measurement on the actual details 
of the multipole structure of the radiative transition. Problematic points in previous calculations of 
the Xc decay angular distributions are identified and clarified. 
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I. INTRODUCTION 

The existing and T polarization measurements 
make no distinction between directly produced states and 
those resulting from the decay of higher-mass states. 3/ip 
and T mesons coming from x decays have, in principle, 
very different polarizations with respect to the directly 
produced ones. In fact, directly produced x and directly 
produced J/?/' or T have different angular momentum 
and parity properties, and originate from different par- 
tonic processes. Moreover, the angular momentum com- 
position of the indirectly produced states is influenced by 
the presence of the accompanying decay photon. There- 
fore, Xc and Xb polarization measurements, together with 
the knowledge of how these states transmit their polar- 
izations when they decay, are essential in the understand- 
ing of the observed J/ip and T polarization patterns. An 
improved account of feed-down effects in quarkonium po- 
larization measurements, and calculations, can shed new 
light in the interpretation of the significant discrepancy 
existing today between the theory predictions and the 
experimental data [1 . 

In this paper we examine how the polarization is trans- 
mitted in the decays from P to S quarkonium states. We 
study the angular distributions of the successive decays 
XciXb) J/V'(T)7 and J/i/'(T) ^ £+1-. We discuss 
the sensitivity of these observable distributions to the an- 
gular momentum composition ( "polarization" ) of the de- 
caying X meson and their additional dependence on the 
orbital angular momentum of the photon. As a result 
of the study, we propose a convenient way of measuring 
X polarizations in high-energy experiments, essentially 
independent of the details regarding the photon detec- 
tion and of the magnitude of the higher-order multipoles 
of the radiative transition. This method is valid irre- 
spectively of the production process (hadroproduction, 
photoproduction, etc.). The paper finishes with a criti- 
cal review of previous calculations of the Xc decay angu- 



lar distributions, identifying and clarifying the causes of 
their seemingly contradictory results. 

II. RADIATIVE DECAY AMPLITUDES 

Throughout this paper, we generically denote by V 
the charmonium and bottomonium ^Si states, J/ip and 
T, and by x the ^Pj states, Xcj and Xbj, with j — 1,2. 
Without loss of generality for the discussions in this pa- 
per, we assume that the ^Pj state Xj is produced in a 
single "subprocess" as a given superposition of 3^ eigen- 
states {z being the quantization axis chosen for the x 
angular momentum), 

j 

ni=-j 

with J^lx; j,m) = j{j -\- 1) Ix; j,m) and J^|x; j,m) = 
m Ix; j,m). Notations for axes and angles are shown in 
Fig. [l]^a) . The total angular momentum carried by the 
photon can have any (non- vanishing) value, while its pro- 
jection along the momentum direction of the 7 (and of 
V), the z' axis, can only be k' = +1 or —1, because 
the orbital component has, by definition, zero projection 
along this direction. In other words, the photon angular 
momentum state is an eigenstate of Jz' but, in general, 
neither of nor of 3z ■ It can be represented as a com- 
plete expansion over eigenstates of and Jz, as 

°° hi -X- 1 ' 

1^;^') = E^^^ E^"'(0'*)l7;^-fc). (2) 

where J2|7; Z, k) = l{l + 1) I7; I, fc), Jz|7; l,k)^k I7; I, k) 
and the coefficients (8,$) [2| are the matrix ele- 
ments of the rotation corresponding to the change of 
quantization axis from z' ("natural" quantization axis 
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of the photon) to z (x quantization axis adopted in the 
measurement). 

The amphtude of the radiative transition from the Xj 
state to the V state plus a photon having spin projection 
k' along z' is 

j oo 
m——j l—l 

I 

X ^ 2?iV(0,<i')(V'7fe'; i,m-k,i,k\n\x;j,m). 

k=-l 

(3) 

The matrix element of the elementary transition can be 
parametrized as 

{V-ik' ; l,m - k,l,k\'H\x; j,m) 

(4) 

= H,i {l,m-k,l,k\j,m), 

where we have factored out the /c'-dependent sign, deter- 
mined by imposing that the photon distribution 

WS-(e,$) = J2 \A{xi^Vjk')\^ (5) 

k'=±l 

is parity invariant and using the property 

PL, (tt - e, TT + $) = (- 1)'-^' e^''^'* Vi (9, $) . (6) 

The sums in Eq. [3] only include terms in which the 
Clebsch-Gordan coefficient {l,m — k, I, k \ j, m) is well de- 
fined, i.e. when 

1 < ^ < j + 1 , |m - fc| < 1 . (7) 

The partial amplitudes Hji denote Xj ~^ ^ transitions 
with the emission of a photon of total angular momentum 
I. In the spectroscopic language they represent electric 



and magnetic 2'-pole radiations (dipole, quadrupole, oc- 
tupole, etc.), indicated with E/ and M^, respectively. The 
two types of transitions differ in their parity properties: 
the electric 2'-pole radiation has parity (—1)' while the 
magnetic 2'-pole radiation has parity (—1)'+^. Since the 
'^Pj and "^S*! quarkonium states have opposite parities, 
the only allowed transitions are El (for all x states), M2 
(for xi and X2) and E3 (for the X2)- Hereafter we use 
the short notations 

hi^Hi.i with 1 = 1,2, 

gi^H2,i with 1 = 1,2,3, ^ ^ 

with the normalizations 

hj + hl = gl+gl + gl = 1. (9) 

In short, hi and gi represent, respectively for xi and X2, 
the relative amplitude of the El transition, /i2 and 32 
the corresponding relative amplitudes of the M2 transi- 
tion, and g3 the relative amplitude of the E3 transition 
(only for the X2 case). The hierarchies 33 < 32 < 5i 
and /i2 < hi are expected. In fact, in the generic ex- 
pansion of the radiation field around a system of oscil- 
lating charges in terms of angular momentum eigenfunc- 
tions, the Z-th term vanishes more and more rapidly, at 
large distance from the origin, as I increases. This be- 
haviour reflects the fact that the wavelength of the emit- 
ted photon, = hc/Ey ~ /ic/(0.4GeV) ~ 3 fm, is 
sizeable with respect to the dimensions of the quarko- 
nium (~ 0.4 — 0.7 fm), so that at the typical distance 
scale r — the electromagnetic field is already only 
weakly sensitive to the internal charge and current distri- 
butions of the radiating object. With respect to the first 
non- vanishing term, higher multipole terms, produced by 
more complex charge/current configurations, are there- 
fore foreseen to be increasingly smaller, even if not neces- 
sarily as suppressed as in nuclear 7-ray transitions, where 
the emitted radiation has a wavelength several orders 
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h2 [%] 


92 




<73 [%1 


Crystal Ball g] 


"•^-2.0 


-33: 


hll 

-30 




E760 [5] 




-14: 


t6 




E835 [B] 


0.2 ±3.2 


-9.3: 


H3.9 
-4.1 


2.0lt« 


CLEO T 


-6.26 ±0.67 


-9.3 ± 1.6 


1.7 ± 1.4 



TABLE I. Higher-order photon multipoles in Xc(i,2) J/V'7 
decays. 



of magnitude larger than the nuclear dimensions. The 
study of the cc radiation multipoles addresses aspects of 
the quark model, including the properties of the bound- 
state wave functions and the electro-magnetic properties 
of the charm quark. For example, the relative contri- 
bution of the M2 amplitudes is significantly dependent 
on the corrections to the charm quark magnetic moment 
fic — '^2m~ 0- "^^^ existing /i2, 32 and 53 measurements, 
for the Xc^ are shown in Table [l] The M2 amplitude con- 
tribution is of order 10% for the Xc2 and even smaller 
for the Xci , although the two most precise Xci results are 
incompatible with each other. No experimental informa- 
tion exists for Xb decays. We will discuss in Sect. [V| the 
effects induced by the higher-order multipole contribu- 
tions on the observable angular distributions. 



III. PHOTON DISTRIBUTION 

The angular distribution of the photon direction in the 
X rest frame, as a function of the x angular momentum 
composition {6„i} and of the photon multipole ampli- 
tudes, is obtained by expanding Eq. [5] The resulting xo 
distribution is spherically symmetric, reflecting the rota- 
tional invariance of the j = angular momentum state 
and the imposed parity invariance of the decay. As for 
the xi decay, the expression of the angular distribution 
is 



w^i(e,$) 



(1 + Ae cos^e 



4^(3 ±Ae) 
+ A$ sin^6cos2$ ± Ae* sin 26 cos$ 
+ A^ sin^esin2$ + Aq^ sin 28 sin$) , 



(10) 



where 



Ae = -(l-3A)[2|6op-(|6+iP±|?>-in], 
A* ---|(l-3A)Re(&;i&_i), 

Ae* = -^(1 - 3A) Rc[6S(6+i - 6_i)] , (11) 
Ai = --|(l-3A)Im(&;i6_i), 



A(!,<, = ^(l-3A)Im[6S(6+i+fc-i)], 



with 



D = 2{1 + A) \bo\' ± (3 - A)(|6+i|' ± 
A = -2/ii /i2 . 



(12) 



The angular distribution of the photon from the X2 
decay, significantly more complex, is 



M^2(e,$) = 



15 



4^(15 ± 5A[^' ±3A[3'^) 



(1 ± A^'^ cos^e 



± A^^ cos^e ± A^'^ sin2ecos2$ ± A^^ sin^ecos2$ 



+ A^^^ sin^ecos4$ + A^'^" sin^esin2$ 



+ A^^^^ sin'*esin2$ ± A^"'" sin*esin4$ 



± A^^^ sin 29 cos $ + X'-^l sin^O sin 29 cos $ 



,(3)_L 
\* 

,(2) 



± X^qI sin'^esin2ecos3<I> ± A^'^-^ sin 29 sin $ 



A^^^-^ sin2esin2esin$ ± A^X sin^esin2esin3$) , 

(13) 



(3)_L 



where 
^(1) 



+ (l-|Ai-|A2)(|6+ip±|6_in 
-(2±iAi-§A2)(|6+2p + |fe-2n], 

Ag^ = f [6I&0P - 4(|5+ip ± + |fe+2p ± \b-2\' 

X'i'> = ^ Re[V6(l ± Ai) b*{b+2 + &-2) 



(2) 
$ 

(3) 



A, 



A. 
A. 

(l)-L 



(3-2Ai -5A2)6;i6_i], 
-ARe[%/6 6S(6+2 ±fe-2) - 4 6;i6_i] , 



^ - ARe(6;2fe-2), 



Im[-%/6(l± Ai)6*(6 



±2 



* D 

+ (3-2Ai-5A2)6;i6_i], 

A^'^^ - Alm[V6b*o{b+2 - 6-2) ± 4 6;i6_i] , 

Ai')^ = AIm(6;26_2), 

^ei - ^ Rc[\/6(1 - |Ai - IA2) bUb+i - 6-1) 

+ (6 ± |Ai - f A2) (6;26+i - bUb-i)] , 
Agi = A Re[V6 55(5+1 - 5_i) - (5;25+i - 5l25_i)] , 
x'il = ARe{b*_,,b_,-bUb+i), 
Ag^-^ = 1 Im[-x/6(l - |Ai - IA2) 5S(5+i ± 6_i) 
+ (6 ± |Ai - f A2) (5;25+i + bUb-i)] , 

-AIm[V6 5o(5+i ± 5_i) ± S^afe+i + 5I2&-1] > 



j,(2)-L 

Agi^ = AIm(5;25_i± 5^25+1), 



(14) 
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FIG. 2. Dependence of the parameters Ae (a) and (b) of the xi photon and dilepton distributions and of the parameters Aq' 
(c) and A^ (d) of the X2 photon and dilepton distributions on the relative contribution of the magnetic quadrupole transition 
amplitude. 



with 



IV. LEPTON DISTRIBUTION 



I? = (10 + Ai - A2)|feol' + (9 - A2)(|fo+iP 

+ (6-iAi + |A2)(|&+2p + |6-2p), 
Ai = 4^2 + 6\/5 5i253 - 2\/55i.g2 - 2gl + Ugig^ 
A2 = 4.9^ + 4^5 3253 + 2%/5.gig2 + Sgl + 4^153 , 

A = 5/(3/?) (A1+A2). 



(15) 



As shown in Fig. [2j the dependence of the photon dis- 
tribution on the X angular momentum configuration is 
very sensitive to the contribution of the higher photon 
multipoles. Figure[2|^c) shows, in particular, that the po- 
lar anisotropy parameter Xq \ at the average value of 52 
measured for the Xc2 (assuming 53 =0), is 30% higher 
than the value expected in the El-dominance case if the 
X2 polarization state is m = ±1 or 70% higher if m = ±2. 
This shows that the derivation of the average polariza- 
tion state in which the x is produced from the observed 
photon angular distribution relies crucially on the knowl- 
edge of the multipole amplitudes. Seen from the opposite 
perspective, we see that the so-called El approximation 
(/i2 = g2 — 93 — 0) is clearly not applicable in the cal- 
culation of the X — ^ decay kinematics expected for a 
given X production mechanism. 



In the parity-conserving case here considered, the gen- 
eral expression for the angular distribution of the dilep- 
ton decay of a vector state is [5] 



w{'&,ip) 



47r(3-t- A^) 



(1 Atf cos^-d 



+ Xip sin 1? cos 2ip + sin 2i9 cos ip 
+ X'^ sin^i? sin 2ip + A^^^ sin sin ip) , 



(16) 



analogous in form to Eq. [TO] The traditional choice 
of axes, adopted in calculations [3 [10] and measure- 
ments [4-7] of the full decay angular distribution for Xc 
mesons produced at low laboratory momentum, is repre- 
sented in Fig. [l]^b), where the V polarization axis, z', is 
the V direction in the x rest frame. With respect to this 
system of axes, any measurement will always find, for in- 
stance in the case of the polar anisotropics for xi and X2 
dileptons (neglecting, for simplicity, the E3 contribution 
in the latter case), the following values: 
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The dilepton distribution in the coordinate sys- 

tem is independent on the % polarization state. This 
choice of axes, while suitable for measuring the contribu- 
tion of the higher-order multipoles, does not provide any 
information on the polarization of the x and, hence, on 
its production mechanism. 

We propose here an alternative definition of the V po- 
larization frame, enabling the determination of the x po- 
larization in high-momentum experiments without the 
need of measuring the full photon-dilepton kinematic cor- 
relations. This definition, shown in Fig. [ijc), "clones" 
the X polarization frame, defined in the x rest frame, 
into the V rest frame, taking the x",y",z" axes to be 
parallel to the x, y, z axes. 

The coefficients of the dilepton distribution can be 
written as a function of the angular momentum com- 
position of the decaying vector state [5], \V) = 

E"=-}an \V; as 



Sm=^i Ix; li'^); the Coefficients of the dilepton an- 

[2|5op-(|6+ip + |&_ip)], 



AA-3|qoP _ 2Re[a;ia-i] 

V2Rc[Qg(a+i - a_i)] 

2 Im[a!j_]^a_i] 

-%/2 Im[ao(a+i -|- a_i)] 



(18) 



where TV = |aoP -|- |a+ip -I- |a_ip. The partial amplitude 
of the Xj decay into a vector state with angular momen- 
tum projection n (= —1, 0, 1) along z" and a photon with 
angular momentum projection k' {— —1, 1) along z' , from 
Eqs.[3]and|4) is 

m=-j 1=1 k=-l (19) 

X (e, $) (-1)^' Hji (1, m - fc, Z, fc I J, m) . 

Inserting these amplitudes into the expressions of the co- 
efficients in Eq. [18] and averaging over the photon states 
/c' = ±1 according to the sum rule 



^(fc' = + l) jj^(fc' = + l) J^{k' = -l) ^{k' = -l) 



X 



3+\ 



Cfc'=+i) 



3+A 



(k' = -i) 



Ar(fc'=+i) 



^(fc' = -i) 



(20) 



3+A 



(fc' = -i) 



with X = Xi,, A^, etc. and iV('='=-i) = iV('='=+i) for par- 
ity conservation, it is possible to obtain the expression of 
the full angular distribution W{Q,^,-d,ip) of the decay 
process X ~^ ^7 ~^ V£^£~. In the following discussion, 
however, we only consider the dilepton distribution, ob- 
tained by integrating 1^(9, ip) over 8 and $. 

In the frame defined in Fig. [l|c), the dilepton de- 
cay distribution of V mesons originating from xa de- 
cays is isotropic. In what concerns the state \xi) = 



gular distribution are: 
1 



A^ 
A<^ 
A^0 



2 

-— Rc(6;i6_i), 
^1 

Di 
2 

Di 



Re[&S(6+i -fe-i)] 
Im(&;i6_i) , 

Im[6S(&+i +&-i)]. 



(21) 



with 



Z?i = D/(l - 35) , 

= 2(1 + 5) |6o|' + (3 - mKi? + I&-1I') , (22) 



5 = i/^r 



The corresponding coefficients for the decay of the 
2' 



state 1x2) = Em=-2 ^™ Ix; 2,m) are 



Aij 

A;p 

A^ip 



with 



D2 
2 

1 

2 

1 



[2|6o|' + |6+i|' + |?>- 



2(|6+2p + |fo-2n] 



Re[\/6 65(6+2 + 6_2) + 36;i5_i] , 

Re[v^6S(6+i - 6_i) + 6 (6;26+i - 6l26_i)] , 

Imh\/6 6*(6+2-6_2) + 36;i6_i], 

Im[-\/6 65(6+1 + 6_i) + 6 (6;26+i + bUb-i)] , 

(23) 



D2^D/{l-5), 

= 2(5 - 5)|6op + (9 - <5)(|6+i|2 + |6_in 

+ 2(3 + 5)(|6+2|2 + |6_2n, 
(5 -2.92 + 152. 



(24) 



Without experimental separation between the xi 
X2 signals (the xo — 1^7 contribution is negligible), the 
dilepton distribution measurement effectively yields the 
corresponding average polarization parameters, implic- 
itly weighted by A^(J=i) and N'^^^^\ respectively the 
numbers of reconstructed dileptons coming from xi and 
X2 decays {X = X^, X^, etc.): 



X = 



3+A 



(3 = 1) 



Af(j = i) 

3+Ay=i' 



jV(3=2) 

3+Ar'' 



(25) 
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V. MEASUREMENT OF x POLARIZATION AT 
HIGH MOMENTUM 

The formulas obtained in the previous two sections 
suggest two remarks. First, with the choice of the 
x"^y",z" axes, the dilepton distribution contains as 
much information as the photon distribution regarding 
the X polarization state. The two distributions are even 
identical when higher-order multipoles are neglected, as 



can be recognized by comparing Eq. [21] with Eq. 11 and 
Eq. [23] with Eq. [l|] for h2 = 92 = 93 0. In this 
limit, for example, = Ae = ^1/3 and +1, respec- 
tively for pure \j,m) = |1,±1) and |1,0) x states, and 
A,9 = Aq' = +1,-1/3 and —3/5, respectively for pure 
|2,±2), |2, ±1) and |2,0) states, while the additional 

(2) (2) 

terms of the photon distribution in the X2 case ( Aq , A^ , 

Aq^ and Aq^) vanish. Second, the dependence of 
the dilepton distribution on the higher photon multipoles 
is negligible, as shown in Fig. [2jb,d) for A^. 

The definition of the x, y, z axes (and, therefore, of the 
x" , y" , z" axes) uses the momenta of the colliding hadrons 
as seen in the x Test frame, so that it requires, in general, 
the knowledge of the photon momentum. However, for 
sufficiently high (total) momentum of the dilepton, the x 
and V rest frames coincide and the x" , y" , z" axes can be 
approximately defined using only momenta seen in the V 
rest frame. For example, if the x polarization axis (z) is 
defined along the bisector of the beam momenta in the 
X rest frame (CoUins-Soper frame |llj). the correspond- 
ing z" axis is approximated by the bisector of the beam 
momenta in the i/ip / T rest frame. The relative error 
induced by this approximation on the polar anisotropy 
parameter is 



(26) 



where AM is the x — mass difference and p is the to- 
tal laboratory momentum of the dilepton. Therefore, for 
not-too-small momentum the frame definition we pro- 
pose coincides with the frame defined in the measure- 
ment of the polarization of inclusively produced i /ip / T 
mesons (Colhns-Soper or helicity, for example). In other 
words, the measurement of the dilepton distribution at 
sufficiently high laboratory momentum provides a direct 
determination of the x polarization along the chosen po- 
larization axis. This determination is cleaner than the 
one using the photon distribution in the x I'cst frame, 
because it is independent of the knowledge of the higher- 
order photon multipoles. 

The above-mentioned approximation, in which the sys- 
tem of axes x" , y" , z" is set without any knowledge of the 
photon momentum, becomes rapidly invalid as p -t> 0. 
The xo case, although of little practical importance in 
the scope of this paper (the branching ratio of the Xco 
decay to i/ij^, for example, is only ~ 1%), can be used 
to give a simple illustration of what happens going from 



AAfl 


= 






A^ 







high to low p. As discussed above, at high p the polariza- 
tion of dileptons from xo vanishes in the helicity frame 
(as well as in any other frame defined ignoring the photon 
momentum), mirroring the perfect isotropy of the pho- 
ton emission in the xo I'cst frame. On the other hand, 
the IS" state coming from the j — Q x state has an in- 
trinsic spin alignment, always opposite to the one of the 
photon (Jv -f = J^j, = 0); in other words, a fully 
transverse 15 polarization is observed if the direction of 
the photon in the xo rest frame [z' in Fig. [ijb)] is taken 
as reference axis. In the low-p limit, when the xo tends 
to be produced at rest in the laboratory, that direction 
tends to coincide with the center-of-mass helicity axis. In 
short, if we choose the center-of-mass helicity frame, the 
V polarization equals the zero polarization of the xo only 
at high momentum, while it changes to fully transverse 
at low momentum, where it simply reflects the intrinsic 
photon polarization. This example shows that the pos- 
sibility to measure the x polarization from the dilepton 
distribution ignoring the photon momentum is strictly 
limited to a kinematic domain where p 3> AAf . How- 
ever, the error in Eq. 26 is already as small as 1% when 



p > A GeV/c, a condition fulfilled, in particular, by es- 
sentially all the quarkonium events collected by the LHC 
experiments. 

The parameters of the dilepton distribution at high 
momentum (Eqs. 21 and 23 with 5 — Q) satisfy charac- 
teristic inequalities. In the xi case. 



-^<A.<+1, |A,|<i^ 

^(a^-^) +6A2^<1, 

,^ , V3/, 1 

\Mip\ < ^ I A^p + - 

(6Ay - if + 6A^^ < 1 for A^ > 



(27) 



In the X2 case, 
5 



16 V^"" 5 



+ 



(28) 



These inequalities continue to be valid in the presence of 
a superposition of production processes leading to differ- 
ent angular momentum compositions of the Xj (see the 
analogous discussion in Ref. [8^ for the direct production 
of a vector state). The corresponding parameter domains 
are represented in Fig.[3| compared with the most general 
constraints valid for vector states, directly or indirectly 
produced. 



VI. COMMENT ON PREVIOUS 
CALCULATIONS 

The angular distributions of the cascade decays Xc — ^ 
3/ip^ were calculated in Ref. [S] (OS) and in 
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FIG. 3. Allowed regions for the angular parameters of the 
dilepton distributions produced by the decay of vector states 
of any origin (light-shaded [S]), of X2 daughters (darker) and 
of xi daughters (darkest). 



Ref. [To] (RSG) for the specific case of low-energy pp 
collisions, where, due to helicity conservation, the Xc is 
only produced in pure eigenstates with eigenvalues 
TO = ±1 (xci) or ±1,0 (xc2)- The two calculations use 
the J/t/i momentum in the Xc rest frame as quantization 
axis for the dilepton, as in Fig. [ijb) , and provide the full 
angular distribution of the correlated photon and lepton 
directions. The result of RSG contradicts the one of OS, 
pointing to a seemingly wrong sign in the last terms of 
the Xc2 distribution (Eq. 10 of OS, corrected into Eq. 20 
of RSG) and of the xa distribution (Eq. 15 of OS, Eq. 27 
of RSG). 

We checked these calculations in two ways, by repeat- 
ing the steps described in the two papers and by com- 
paring them to our own calculation for the full decay 
distribution in the special case of pure iz eigenstates. In 
the latter case, we have applied a rotation of the lepton 
variables from the x" ,y" ^z" system adopted in our cal- 
culation to the a:', y', z' system adopted in OS and RSG. 
We found that, except for an apparent misprint of OS 
(the fifth line of Eq. 11 in OS has a wrong numerical co- 
efficient, corrected in Eq. 21 of RSG), both calculations 
are correct. RSG argued that OS used two inconsistent 
conventions for the reduced rotation matrices d\p adopt- 
ing one ordering of the indices i and ] (the one used in 
RSG) in the description of the J/-!/) — >■ l^t~ process and 
the reverse ordering in the description of the Xc — J/V" 7 
process. We have verified that, instead, the conventions 
are everywhere consistently used, while RSG did not con- 
form to the calculation of OS and adopted a different def- 
inition of the photon angle. OS refers, for the adopted 



notation, to Ref. [T^, where the axes definitions are de- 
scribed in the first figure of the paper. Even if there is 
no explicit mention in the text, the angle d in the figure 
(which we denote by O in our Fig. 1) is, unmistakably, 
the angle formed by the photon momentum with the an- 
tiproton direction in the Xc rest frame, while d' (which 
we denote by in our Fig. 1) is the angle formed by 
the lepton momentum in the J/'0 rest frame with respect 
to the J/'f/' momentum in the Xc rest frame. RSG uses 
the same definition of 0', but an opposite definition of 
B: "We will work in the xj rest frame with the Z axis 
taken to be in the direction of ?/;. The p direction is in 
the X-Z plane, making an angle Q with the Z axis" . As 
a consequence, when a certain reduced d matrix is used 
in OS to rotate the quantization axis by an angle 0, the 
inverse rotation must appear in the calculation of RSG. 
If dlj{9) represents a given rotation, the inverse rotation 
can be denoted either by exchanging i with j (this in- 
duced RSG's misinterpretation of the discrepancy) or by 
replacing with 27r — 9. This explains the different sign 
in the term proportional to sin 29 resulting from the two 
calculations. The remaining terms, depending on cos^0, 
are not sensitive to such a redefinition of the angle. 

In short, each of the two calculations is correct, if they 
are made with the matching angle definition. If, on the 
contrary, the definition of 9 used by OS is used together 
with the distributions functions derived in RSG, or vice- 
versa, a wrong sign appears in the term proportional to 
sin 29, leading to unphysical results. In fact, this artificial 
change of sign is not reabsorbed in a different definition of 
sign and/or magnitude of the higher-order multipole am- 
plitudes: already in the El approximation, the physical 
correlation between photon and lepton angles is substan- 
tially altered by such a mistake. To evaluate the impor- 
tance of this problem, we assumed the angle definitions 
of OS and used the formulas derived in RSG, transposing 
them, by rotation, to the system of axes used in our cal- 
culations [Fig. [ijc)]. As a result of this forced mistake, 
we arrive to a physical result which is almost opposite to 
the correct one: the lepton distribution, instead of be- 
ing a perfect clone of the photon distribution (in the El 
limit), becomes a consistently smeared, almost isotropic 
distribution, for whatever polarization state of the x (in 
other words, the domains of the xi and X2 dilepton pa- 
rameters, represented in Fig.jSj are reduced to small areas 
around the origin). 

We have noticed that the measurements of E760 [5] 
and E835 |6], included in the present world averages of 
/i2 , 52 a-nd in the Review of Particle Physics '13] , seem 
to be affected by this kind of misunderstanding. Both 
analyses define the photon angle 9 as "the polar angle 
of the J/ip with respect to the antiproton", as in OS, 
but the formulas are taken from RSG (Table II in the 
E760 paper and Tables IV-V in the E835 paper repro- 
duce Eqs. 20 and 27 of RSG). On the other hand, the 
quality of the global fits of the data using the adopted 
parameterization is rather good and the measurements 
of the higher-order multipoles are compatible with the 
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CLEO results f7|, suggesting that the inconsistency be- 
tween formulas and angle definitions might simply be an 
editing mistake in both experimental papers. 



VII. CONCLUSIONS 

We have derived the expressions of the angular dis- 
tributions of the radiative decay from a ^Pj state to a 
^Si state and of the dilepton decay of the latter. No 
selection rules specific to certain quarkonium production 
mechanisms have been used and the choice of the polar- 
ization frame for the directly produced states has been 
kept completely general. 

We have shown that the x polarizations can be mea- 
sured (for not-too-low-momentum experiments) directly 
from the angular distribution of the dilepton decay in 
the i/ij) / T rest frame, with respect to the same kind of 
system of axes (Collins-Soper, helicity, etc.) adopted in 
inclusive J/^ / T measurements. 

In fact, the dilepton distribution in the i/ip / T rest 
frame is a clone of the photon distribution in the x rest 
frame, stripped of the contribution of the higher multi- 
poles of photon radiation. 

This represents a significant advantage, given that such 
contributions — measured to be quite important in the 
Xc2 case, poorly known (due to contradictory measure- 
ments) in the Xci case, and still unmeasured for the bot- 
tomonium family — can have a very large impact in the 
measurement. Furthermore, a simultaneous determina- 
tion of X polarization and of the multipole parameters is 
scarcely feasible at hadron colliders. 

An additional advantage of this method is that it 
does not use the photon measurement to reconstruct the 



event-by-event decay topology. This means that, con- 
trary to previous expectations, the measurement of x po- 
larization is not intrinsically more challenging than, for 
instance, the measurement of the Xci/Xc2 cross section 
ratio (a measurement presently being done by several ex- 
periments at the Large Hadron Collider). In both cases 
the analysis needs to identify an event sample where the 
J/'0 (or T) dilepton is associated to a photon giving an 
invariant mass of the fi^fi^j system in the x mass re- 
gion. This is usually done using photons reconstructed 
by the conversion method, given that the tracking of the 
electron-positron pair gives good enough resolutions to 
resolve the Xci and Xc2 resonances. Naturally, a larger 
event sample is needed for a multi-dimensional angu- 
lar analysis. But there are no extra difficulties related 
to photon backgrounds or reconstruction efficiencies de- 
pending on the decay angles, specific to the measurement 
of the polarization (as would be the case using the pre- 
viously available methods). 

It is worth reminding that a certain numerical value 
of the observable polarization parameters corresponds to 
very different quantum-mechanical states of the xi and 
X2 (e.g., — +1 can reflect the Jz = state of the xi 
or the Jz = ±2 state of the X2)- Therefore, a reliable ex- 
perimental discrimination between the 3/tp or T coming 
from the decays of these two states is crucial for a proper 
understanding of x polarization. 

We have also pointed out misunderstandings in pre- 
vious calculations, which may have affected some of the 
existing measurements of the higher-order photon multi- 
poles in Xc decays. 
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